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Abstract 

We propose a new scenario for early cosmology, where an inflationary de Sitter 
phase is obtained with a ghost condensate. The transition to radiation domi- 
nance is triggered by the ghost itself, without any slow-roll potential. Density 
perturbations are generated by fluctuations around the ghost condensate and 
can be reliably computed in the effective field theory. The fluctuations are scale 
invariant as a consequence of the de Sitter symmetries, however, the size of the 
perturbations are parametrically different from conventional slow-roll inflation, 
and the inflation happens at far lower energy scales. The model makes definite 
predictions that distinguish it from standard inflation, and can be sharply ex- 
cluded or confirmed by experiments in the near future. The tilt in the scalar 
spectrum is predicted to vanish (n s = 1), and the gravity wave signal is negligi- 
ble. The non-Gaussianities in the spectrum are predicted to be observable: the 
3-point function is determined up to an overall 0(1) constant, and its magnitude 
is much bigger than in conventional inflation, with an equivalent /nl — 100, not 
far from the present WMAP bounds. 



1 Introduction 



Inflation is an attractive paradigm for early cosmology. An early de Sitter phase of the 
universe elegantly solves the horizon and flatness problems. Further imagining that this 
phase is driven by a slowly varying scalar field with a flat potential gives a beautiful 
mechanism for generating the observed density perturbations, which have their origin in 
the quantum fluctuations of this light field during inflation. This minimal picture of slow- 
roll inflation predicts the fluctuations to be nearly scale invariant and Gaussian, in excellent 
agreement with observation. 

While this qualitative picture is very attractive, it is frustratingly difficult to make much 
sharper predictions that could definitively exclude or confirm it in future experiments. 
There are nevertheless some generic expectations for what should be seen. For instance, 
in most models, the tilt in the scalar spectrum is non- vanishing, and is typically of order 
n s — 1 ~ (1/N e ) ~ 2 x 10 -2 , where N e is the number of e-folds to the end of inflation when 
typical cosmological scales leave the horizon. The current bounds on n s — 1 are getting 
close to this size, while the Planck satellite will improve the sensitivity to ~ 3 • 10 -3 . 
Slow-roll inflation leads to the expectation that Planck should find n s ^ 1, but this is 
not a hard prediction, and if Planck does not find a deviation from n s = 1 we can not 
definitely exclude it. If a tilt is found, it would be interesting to look for finer structures in 
the power spectrum, such as bumps or wiggles of various kinds. However, one can always 
engineer bumps and wiggles in the inflaton potential to reproduce these features in the 
power spectrum, so there is nothing here either that can exclude or confirm the standard 
picture with certainty. The observation of a gravitational wave signal in the CMBR would 
be extremely exciting and it would indicate that the energy density driving the accelerated 
expansion is very high, of order Mq UT . Minimal slow roll inflation predicts a relationship 
between the tilt of gravity wave spectrum n g — 1 and the scalar one n s — 1, and if this is 
found it would be a confirmation of the standard picture. However, we would have to be 
lucky for the gravity wave spectrum to be observable in near future experiments, to say 
nothing of the tilt in the gravity wave spectrum, and if the inflationary prediction is found 
not to hold, it would by no means exclude the theory; deviations from this relationship can 
come about by integrating out new physics between the Planck and Hubble scales during 
inflation [TJ[2]. 

The minimal picture of slow-roll inflation does make one very sharp prediction, how- 
ever: the spectrum is predicted to be Gaussian to a very high degree. For single field 
inflation models, the 3-point function was recently calculated to leading order in the slow- 
roll parameters by Maldacena The dimensionless size of the effect was found to be 
parametrically at least of the same order as Sp/p itself, (and in fact further suppressed by 
slow-roll parameters). The prediction is very firm, since the leading non-Gaussianities are 
generated by the interaction of the inflaton with gravity, and are therefore determined by 
the symmetries of the theory. However, the magnitude is about a factor of ~ 100 smaller 
than what might be observable in the future experiments. Because of this firm prediction, 
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however, non-Gaussianities are the best and most exciting place to look for deviations from 
the standard picture of slow-roll inflation. Indeed there are by now a number of different 
ways of generating density perturbations from sources other than the fluctuations of the 
slowly rolling inflaton, and in these models it is possible to get larger (and observable) 
non-Gaussianities [UG3E]- Any observation of non-Gaussianities necessarily implies the 
presence of additional light fields and a non-standard mechanism for generating density 
perturbations. 

Given the difficulty of conclusively verifying the standard inflationary paradigm, it is 
worth looking for alternatives, especially if they are more predictive in some way. Given 
the ease with which a de Sitter phase of expansion obliterates the horizon and flatness 
problems, it would be nice to keep this, but is there an alternative to the standard slow- 
roll picture? 

In this paper we will describe such a model. Recently, a new possibility of having 
de Sitter phases in the universe in a way differing from a cosmological constant has been 
proposed which also modifies gravity in the infrared in a non-trivial and consistent 
way. It can be thought of as arising from a derivatively coupled ghost scalar field (f> which 
"condenses" in a background where it has a non-zero velocity 

(</,) = M 2 -> (<f>) = M 2 t . (1.1) 

The novelty here is that, unlike other scalar fields, the velocity <fi does not redshift to 
zero as the universe expands, it stays constant, and indeed the energy momentum tensor 
is identical to that of a cosmological constant. However, the ghost condensate is not a 
cosmological constant, it is a physical fluid with a physical fluctuations ir defined as 

(f) = M 2 t + iT. (1.2) 

The ghost condensate then gives an alternative way of realizing de Sitter phases in the 
universe. The symmetries of the theory allow us to construct a systematic and reliable 
effective Lagrangian for ir and gravity at energies lower than the ghost cut-off M. Neglect- 
ing the interactions with gravity, the effective Lagrangian for n (around flat space) has the 
form 

where a and (3 are order one coefficients. Note that this Lagrangian is non-Lorentz in- 
variant, as should be expected, since the background <fi = M 2 breaks Lorentz invariance 
spontaneously (the it field can be thought as the Goldstone boson for this symmetry break- 
ing). The low-energy dispersion relation for ir is of the form 

- 2 = ° 2 ^- ("> 

It is straightforward to couple this healthy sector to gravity, which leads to a variety of 
interesting modifications of gravity in the IR, including antigravity and oscillatory modu- 
lation of the Newtonian potential at late times and large distances jjj. 



2 



Given that this theory provides an alternative way of having a de Sitter phase, it is 
natural to try to use the field eft as the inflaton. Since we have a different value for the 
Hubble scale during inflation than the present, the shift symmetry on <p must be broken. 
However, it is technically natural to assume that there is a potential V(4>) which has the 
form e.g. V((f>) = Vq for (j) < 0, V{4>) = for (p > 0. The shift symmetry is still good for 
4> < and 4> > 0, so that this shape is radiatively stable. Of course this is an exaggeration 
of a smooth shape for the potential, and we could have variations on each side over scales 
smaller than the cutoff M of the ghost theory. The theory inflates for <p < 0, and evolves to 
radiation dominance for (j) > 0. Note that we need not assume that Vq itself is of order M 4 . 
For instance, as in hybrid inflation, <f> could couple to matter with shift-symmetry breaking 
couplings around = 0; these couplings can trigger a phase transition at 4> = and Vq 
can be the energy of the false vacuum in the matter sector. We will not concern ourselves 
with describing any specific model for the transition or with justifying the different energy 
scales, M and Vo, of the model; it suffices to say that the set-up is completely technically 
natural, i.e. radiatively stable. From this minimal set-up, however, we can make definite 
and experimentally relevant predictions for the density perturbations, which is our main 
focus in this paper. 

Let us begin with a qualitative picture and estimates for the size and level of non- 
Gaussianity in the perturbations, and later derive all the results more formally. There are 
two important differences here from ordinary inflation. First, there is no slow-roll. Even in 
the approximation where the potential for <fi < is exactly flat, <fi = M 2 is non-zero, and (f> is 
continuously heading towards = where inflation ends, continuing with about the same 
velocity afterwards. As we will see in a moment, density perturbations are generated here 
from the fluctuations of the tt field, and are again scale invariant as in inflation. However, 
since there is no need for any slope in the potential, there is no reason to expect any tilt 
in the scalar spectrum. 

The second important difference with standard inflation concerns the size of the fluctu- 
ations in (or equivalently tt). Since the effective Lagrangian for tt is non-relativistic, and 
in particular there are no k 2 spatial kinetic terms, the fluctuations of tt are less suppressed 
than in a relativistic theory. In a relativistic theory, the fact that scalar fields have scaling 
dimension 1 tells us that the size of the fluctuations of a scalar field inside a region of size 
R is given by ~ 1/-R, similarly at a frequency E it is given by E. In usual inflation, the 
inflaton fluctuations freeze when they have a typical energy E ~ H, so that their typical 
size is 5(j) ~ H. We can determine the size of the fluctuations in our case by a simple scaling 
argument familiar from power-counting for non-relativistic effective theories (HJ- Suppose 
we scale energies by a factor of s, E — > sE, or alternatively t — » s~ 1 t. Clearly, because of 
the ijj 2 oc A; 4 dispersion relation, we have to scale k differently, k — > s 1 / 2 k or x — » s~~ x l 2 x. 
We then determine the scaling dimension of it by requiring the quadratic action to be 
invariant, and we find that tt has scaling dimension 1/4 

tt -> s 1/4 vr . (1.5) 
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Now, -/r has mass dimension 1, so the fluctuations at frequencies of order the cutoff M is 
5ttm ~ M. But the fact that tt has scaling dimension 1/4 tells us that the fluctuation at a 
lower energy E is Site ~ (-EM 3 ) 1 / 4 . In particular, the size of ghost fluctuations that freeze, 
as usual, by Hubble friction when its frequency in of order E ~ H, is 

8tt h ~ (HM 3 ) 1 / 4 . (1.6) 

Of course, for consistency of the effective theory, we must have H <C M; note that as 
expected these fluctuations are much larger than H in this limit. 

Let us now estimate the size of the density perturbations. As usual, <fi fluctuates and is 
stretched out until the Hubble damping becomes important at frequency E ~ H; note that 
this does not correspond to k ~ H but rather, from the dispersion relation, k ~ \J ELM. 
The fluctuation bixn causes inflation to end at slightly different times in different places, 
and so as usual we have the estimate 

S -£~HSt=**™. (1.7) 

P <\> 

Now both of the differences with the usual story come into play. First, <f> = M 2 , having 
nothing to do with slow-roll parameters. Second, Sttjj is much larger. We then find 



which can be compared with the usual inflationary case 

Sp H 
7 ~ M Fl yfe ' 



(1.9) 



where e is the usual slow-roll parameter: e = • {V'/V) 2 . Since d> continues to 

have the same velocity today, leading to the IR modifications of gravity described in [7j, 
M is limited to be at most 10 MeV or so: we must have H ~ 1 keV at largest for 
Sp/p ~ 10 -5 . We therefore generate the correct magnitude density perturbations at far 
lower energy scales than usual inflation. Indeed, the maximum energy scale is of order 
Vq ~ H 2 Mp\ ~ (1000 TeV) 4 . This implies that the gravity wave signal is completely 
negligible. 

We can then reproduce the spectrum of perturbations, with the sharp predictions that 
(a) \n s — 1| = and (b) there are no gravitational waves. The observation of any of these 
would decisively rule out our model. 

Things get more interesting when we look at the non-Gaussianities in the model. Be- 
cause all the scales are so much smaller than the Planck scale, there is no effect coming 
from gravitational interactions. The dominant effect comes from the trilinear interaction 
in the tt effective Lagrangian 

A^(Vvr) 2 . (1.10) 
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This leads to a non-zero three-point function for the density perturbations. To get an idea 
of the dimensionless size of this effect, we need to find the dimensionless size of this coupling 
at an energy of order H. We can do this easily since we know the scaling dimension of 
t, x and 7r: we find that the coefficient of this operator has scaling dimension 1/4; it is 
an irrelevant operator, but just barely! The dimensionless size of this interaction at the 
cutoff M is ~ 1; scaling it down to an energy of order H we get an estimate for the 
non-gaussianity of the perturbations 



These are much larger than in standard inflation, where NG ~ e • (5p/p)\ In fact, this 
estimate suggests that the non-gaussianity may be too large, ~ 10 _1 , but in fact the com- 
plete analysis we do below shows that it is naturally just at the level of the present WMAP 
constraints. The estimate clearly shows the physical origin of the effect, as being due to the 
unusual fluctuations and scaling dimensions of interactions in the ir effective theory. Fur- 
thermore, this operator is the dominant (least irrelevant) operator at low energies. As such, 
the dependence of the three-point function on momenta can be unambiguously predicted 
up to small corrections, as we will see in detail below. 

Before turning to the detailed analysis of our scenario, we point out that the idea of 
an accelerated expansion of the universe driven by a non-minimal kinetic term was already 
discussed under the name of k- inflation j^j . K- inflation is based on some assumed form for 
the higher derivative terms, which can only be justified by the knowledge of the fundamental 
theory. On the contrary our model is based on a consistent effective field theory around 
the ghost condensate, so that all the predictions are not UV sensitive. 



Quantum fluctuations of the ir field in de Sitter space become, as usual, classical fluctu- 
ations at large scales, when they freeze and become constant. In ref. 7_ it was shown 
that the gravitational potential in longitudinal gauge, $ = ^, decays to zero outside the 
horizon. This means that 7r fluctuations do not gravitate at superhorizon scales in the 
pure de Sitter background. This is quite intuitive: in the limit in which the shift symmetry 
7r — > ir + const, is exact, different parts of the Universe with a different value of ir are 
equivalent, so that the metric is not perturbed. 

This is rather different from what happens in standard inflation, where the surfaces 
of constant inflaton are also of constant energy density and the gravitational potential, 
therefore, does not decay to zero outside the horizon. Our model is similar to the usual 
slow-roll inflation in the limit in which the potential becomes completely flat, keeping the 
produced density perturbations constant: e — > with H/{M-p\y/e) = const. In this limit 
the comoving surfaces of constant inflaton have the only meaning of setting the time to 
the end of inflation, but there is no metric perturbation before reheating: — > as e — > 0. 





2 Generation of density perturbations 
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The calculation of density perturbations will therefore proceed as usual: we can calcu- 
late the gauge invariant quantity £ ^Uj, proportional to the spatial curvature of comoving 
surfaces through the position dependent time shift 

c = -5*. (2-i) 

9 

This quantity will remain constant outside the horizon, independently of the details of the 
reheating process and it will seed all the perturbations we observe today. For a simple 
check of the conservation of £ in the limit of instantaneous reheating, see the Appendix. 

The situation is in some sense simpler than in conventional slow-roll inflation. In that 
case we are forced to switch to the £ variable at horizon crossing, because we know that ( is 
constant outside the horizon, while the inflaton perturbation will keep on evolving because 
of the non-zero potential. The evolution of the inflaton perturbation outside the horizon 
is such that it cancels the time evolution of H/<fi in IJ2.1JI to give a constant £. In our case 
both £ and ir are constant outside the horizon before the end of inflation as H and <fi are 
time-independent. So we can also easily follow the surface of constant ir until reheating. 

Let us now calculate the spectrum of density perturbations. In presence of the ghost 
condensate gravity is modified in the IR and this modification is characterized by a typical 
time scale T -1 , with T ~ M 3 /Mp V and a typical length scale to -1 , with m ~ M 2 /Mpi (Jj. 
The two scales are different because the condensate obviously breaks Lorentz invariance. 
To avoid large modification of gravity nowadays we have to require T < Hq, where Hq is 
the present Hubble constant [7] . Obviously this implies that gravity is not modified during 
inflation, i.e. 

r«m«ff. (2.2) 

This relation is equivalent to the decoupling limit Mp\ — ► oo, keeping H fixed. Therefore 
we can study the ir Lagrangian alone, without considering the backreaction on the metric 
through the Einstein equations. At linear level we get 

a 2 fk 4 \ 

TTfc + 3H7T k + f -j J TT k = 0, (2.3) 

where k is the comoving wavevector and a the scale factor. This equation describes a free 
field in de Sitter space, with the modified dispersion relation uJ 1 oc k . 

To calculate the spectrum of 7r fluctuations, we quantize the field as usual 

7T fc (i) = w k (t)a k + w* k (t)al k , (2.4) 

where w k (t) satisfies (|2.3j) . A qualitative study of eq. ()2.3|) tells us that the wavefunc- 
tion goes to a constant when the expansion rate H exceeds the oscillation frequency 
aM~ 1 (k/a) 2 : due to the modified dispersion relation the tt field freezes when its phys- 
ical wavelength is of order (a^HM ) -1 / 2 and not H~ l as usual. 
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We can easily estimate the asymptotic amplitude of Wf-, which gives the spectrum of 
the 7r perturbations. Well before freezing eq. ()2.3j) describes an harmonic oscillator with 
slowly varying frequency and with a friction term. Both the variation of frequency and 
the friction reduce the oscillator energy. Using the equipartition of kinetic and potential 
energy, it is easy to obtain the time variation of the total energy E of the oscillator 

E(t) = —5HE(t) , (2.5) 

which implies that the quantity E{t) ■ a 5 is constant before freezing. We can rewrite this 
using the amplitude of oscillation <$7r(i) and the frequency u>(t) as 

5Tr(t) 2 Lu(t) 2 a 5 = const . (2.6) 

Now we can get the asymptotic amplitude using this conservation law between one time 
well before freezing (where we set a = 1) and the freezing point u(t) ~ H. At the first 
point we use the usual flat-space normalization of the wavefunction Sir ~ cj -1 / 2 and we get 
the final result 

S ^(k) ~ * j^— , (2.7) 

compatible with what we got in the introduction using the scaling dimension of the oper- 
ators. The spectrum of tt fluctuations is scale invariant. This is not a surprise; it is just 
the consequence of time translation invariance in the de Sitter background; every mode 
freezes at a fixed amplitude when it reaches a certain physical size and this is enough to 
give a scale invariant spectrum, independently of the dispersion relation of the tt field and 
the fact that its modes freeze at a scale smaller than the de Sitter horizon. 

The qualitative behavior above can be checked by the explicit solution of eq. (|2.3|) . To 
simplify the equation we can go to conformal time dr\ = dt/a (we take r\ = —(aH)^ 1 during 
inflation) and write the equation for the variable ut = Wk ■ a- We get 

„ fa 2 H 2 k 4 9 2 \ 
<+{- 1 ^^ 2 -- 2 )u k = ^ (2.8) 

where the derivatives are now taken with respect to ry. The solution of eq. (|2.8j) with the 
correct flat space limit for very short wavelength (i.e. r/ — ► — oo) is given by 

u k (v) = ^V^HM(E^ v A , ( 2 . 9 ) 



3 / 4 V 2M 

where H is the Hankel function, linear combination of Bessel functions: Hp {x) = J v {x) + 
iY v {x). The spectrum of tt can be calculated from the asymptotic behavior of (J2.9|) in the 
limit T) — > 0. The Bessel function J_3/4 is dominant and we get 

iV2^ ( M \ 3/4 i , 

"•'("^TWiiUJ k " ■ <2 ' 10) 
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The 7r spectrum is therefore given by 

* = «" - 0)| 2 = ffra(1/ ^ • (2.11) 
From this we get the primordial curvature spectrum through (|2.1j) 

& = » i ^yv3 /4 . (, 12) 



v^r(l/4) \AT(V4) 
The COBE normalization fixes P^ /2 - 4 - 8 • 10 ~ 5 - 

3 Non-gaussianities 

One of the sharpest prediction of inflation, at least if the inflaton itself is responsible for 
the generation of density perturbations, is that these density fluctuations should have a 
probability distribution very close to gaussian. There are two reasons behind this property. 
First of all the inflaton field must have a very shallow potential which allows a slow-roll 
phase; this means that the inflaton self-interactions are very small so that it behaves like a 
free field. Non-linearities coming from General Relativity, even if suppressed by the Planck 
scale, are in fact bigger than those coming from the inflaton potential, but anyway too 
small to be experimentally testable: the fluctuations are gaussian up to a level of 1(T 6 0. 
One can imagine that additional self-interactions for the inflaton are coming from higher 
dimension operators with derivatives: these operators are compatible with slow-roll if they 
respect a shift symmetry on the inflaton (f> — > ^+ const. The contribution of these operators 
cannot be very big anyway: the lowest dimension operator, (V</>) 4 , has scaling dimension 
2, so that its contribution is suppressed by (H/A) 2 , where A is the cut-off scale. As A 
cannot be too small if we want the classical motion of the inflaton to be under control (i.e. 
cf> <C A 2 ), we get a contribution which cannot be much bigger than 10 -5 j!2j . 

In our model we do not have any potential for <j> during inflation and the only non- 
linearities come from higher dimension operators. As discussed in the introduction, these 
interactions are actually much more relevant than in standard inflation, because of the non- 
relativistic scaling of the operators. The most relevant (or better least irrelevant) operator 
is 

V(W) 2 . (3.1) 



M 2 

It has scaling dimension [energy] 1 / 4 , all the others have dimension bigger or equal to 1/2 
so that they will be less important in calculating the non-gaussianities of the density 
perturbations. The very small dimension of this operator tells us that this will be much 
more important than higher dimension operators in standard inflation, which have at least 
dimension 2. Again we implicitly assumed a decoupling limit, neglecting the backreaction 
of 7T on the metric: in this way we neglect all the non-linearities intrinsic to General 
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Relativity, which are suppressed by Mpi and therefore subleading with respect to the tt 
self-interactions. 

Let us estimate the size of the non-gaussianities comparing the non-linear correction 
tt(Vtt) 2 with the linear terms in eq. (|2.3|) . close to the freezing point u ~ H. The corrections 
will be given by the ratio 



h 2 p w - v^r(i/4) \m) 



a -r/4 „ 3 . w -2 . p . a -8/5 5 (3 2) 



where in the last equality we used the COBE normalization eq. (|2.12|) . Parametrically we 
obtain the same estimate we got in the introduction using the scaling dimension of the 
operator. Taking a ~ P ~ 1 the level of non-gaussianity is quite substantial and close to 
what is allowed by experiments. 

Despite the order one uncertainty we expect that the non-gaussian contributions are 
experimentally detectable. The angular distribution of the 3-point function clearly de- 
pends on the wavefunction of the tt field and on the specific trilinear interaction and could 
in principle be used as a smoking gun of the model. We therefore turn to the explicit 
calculation of the 3-point function given by the interaction 7r(V7r) 2 . 

Before doing the explicit calculation we must clarify the relation between tt and the 
gauge invariant quantity £. At the end we are obviously interested in the 3-point function 
of £, a s this quantity remains constant outside the horizon (also at non-linear order 
and it is directly related to the fluctuations we observe today. In doing the calculations it is 
much easier to calculate the 3-point function of the tt field and only at the end turn to the 
£ variable. Physically we calculate the 3-point function of the tt field after its fluctuations 
are frozen outside the horizon then, before the end of inflation, we switch to comoving 
coordinates where tt is constant on surfaces of fixed time, as tt fixes when inflation ends at 
each point. The position dependent time shift between the two sets of coordinates is given 
by (|2.1j) . It is very important to note that this relation is valid also at non-linear order in 
tt as 4> an d H are constant so that they do not depend on tt. This is different from what 
happens in slow-roll inflation where the relation between £ and the inflaton fluctuations is 
non-linear pj . In our case there are no additional non-linearities going from the variable tt 
to the two 3-point functions are just proportional through the relation (|2.1|) . 

In the calculation of the 3-point function we start from the interaction Lagrangian 



Ant 



t" (W)2 ] • (3 - 3) 



We are interested in the value of the (tttttt) correlator at a time t when the fluctuations are 
frozen outside the horizon, so we have to evolve it starting from an initial time deep inside 
the horizon 



{^kA t )^k 2 { t )^k a {t)) = -i d£ ^ 



t 

3„ 



TTfclWTTfcaC*) 71 "^^)' / d X Kindt' 



(3.4) 
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The expression can be evaluated using the decomposition <|2.4f) of the it field, giving 

i 

1 



oc 



dVjj-wl^wl^w^ivKh ■ k 2 ) + symm. + c.c. 



The sum over the two analogous expressions with the derivative acting on the other wave- 
functions is indicated with '+symm'. The integral must evaluated with the prescription 
that the oscillating functions become exponentially decreasing, i.e. below the negative 
real axis: this fixes the vacuum of the interacting theory at rj = — oo. Using the explicit 
expression for the wavefunctions ()2.9I) we get to the final expression 

^ k ^ = -T(TW^\'^) (27r) <5 (^» (3 - 6) 

J drj r]- 1 F*(ri)F* (j^VJ F> * {j^ 71 ^ h(h ■ k 2 ) + symm. + c.c. 

where 

F[x) = f-{-x)^H^/2) . (3.7) 
We can finally translate the result into the £ variable, using eq. (|2.1[) 

= (2^ 3 (E h) (3.8) 

—^ 2 Re y drj ri~ 1 F*(rj)F* (j^vj F '* (j^ 7 ^) ^(fa ' ^ + s y mm - 

As there are no poles one can rotate the contour of integration along the direction 
oc (— 1 — i) so that it converges exponentially. 

In the analysis of the data (see e.g. ^3]) it is usually assumed that the non-gaussianities 
come from a field redefinition 

C(x) = C 9 (x) - lhUC 2 g (x) - <C 9 2 » , (3-9) 

where ( g is gaussian. This pattern of non-gaussianity, which is local in real space, is 
characteristic of models in which the non-linearities develop outside the horizon. This 
happens for all the models in which the fluctuations of an additional light field, different 
from the inflaton, contribute to the curvature perturbations we observe. In this case non- 
linearities come from the evolution of this field outside the horizon and from the conversion 
mechanism which transforms the fluctuations of this field into density perturbations. Both 
these sources of non-linearity give a non-gaussianity of the form because they occur 
outside the horizon. Examples of this general scenario are the curvaton models jH], models 
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with fluctuations in the reheating efficiency [Hj and multi- field inflationary models @|. In 
the data analyses (|3.9j) is taken as a simple ansatz and limits are therefore imposed on the 
scalar variable /nl- The angular dependence of the 3-point function in momentum space 
implied by (|3.9|) is given by 

(C kl Ck 2 Ck 3 ) = (2nf5 3 (^) (2^) 4 ( " pNlPl) • (3.10) 

In our case the angular distribution is much more complicated than in the previous expres- 
sion so the comparison is not straightforward. 

We can nevertheless compare the two distributions (|3.8|) and (|3.1Uj) for an equilateral 
configuration and define in this way an "effective" /nl for k\ = k% = k^. Evaluating 

^ 1 /2 

numerically the integral in (|3.8j) for k\ = ki = k^ and using COBE normalization ~ 
4.8 x 1CT 5 and eq. l|TT2j) we get to 

/X 1 -^85-/3-a- 8 / 5 . (3.11) 

The result is rather suppressed, taking a ~ ~ 1 with respect to the estimate (j3.2[) . which 
would give /nl ^ 10 3 . There are two reasons tor this suppression: first of all n — > 
near freezing and this suppresses the integral in (|3.8|h Second the wavefunction is slightly 
suppressed with respect to its natural size for 77 — s- and this gives an additional small 
suppression. The result must be considered just as an estimate of the effect because order 
one factors are not under control. On the other hand the angular dependence in ()3.8|) 
is well defined and different from all the models generating observable non-gaussianities 
through an additional light field; the angular pattern represents a potential smoking gun 
of the model. 

The present limit on the non-gaussianity parameter from the WMAP collaboration ^3] 

-58 < /nl < 138 at 95% C.L. (3.12) 
is quite close to the estimate (|3.11|) . 

We expect that the limit would change appreciably using for the analysis the angular 
pattern l(3~H|> . The reason is that our pattern is qualitatively quite distinct from the one 
parametrized by /nl- The behavior in the limit in which one of the wavevector goes to 
zero is, for example, rather different from what happens in IjM.lOJ) . In this limit one of the 
fluctuations has a very long wavelength, it exits the horizon and freezes much before the 
other two and acts as a sort of background. Let us take £3 very small and assume that a 
spatial derivative acts on this background in the interaction Lagrangian (|3.3|) . The 2-point 
function (^1^2) depends on the position on the background wave and it is proportional to 
9j7T3 at linear order. This variation of the 2-point function along the tt^ wave is averaged to 
zero in calculating the 3-point function (^^^3), because the spatial average (^39^3) 
vanishes. So we are forced to go to second order and we therefore expect the integral in 
(|3.8|) to go as making the full 3-point function proportional to k^ 1 , after we multiply 
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for the 7T3 spectrum. Note that taking the time derivative in (J3.3|) to act on the background 
gives a subleading contribution as the wavefunction goes to a constant as fc|. On the other 
hand the "local" behavior (|3.10j) goes as fc^ 3 in the limit we have considered. In our model 
the 3-point function is generated by a derivative interaction which favors the correlation 
of modes freezing roughly at the same time, while the correlation is suppressed for modes 
of very different wavelength. The same happens in standard inflation when we study the 
non-gaussianity generated by higher derivative interactions ^2]. The implication of very 
different angular patterns for the experimental limits on non-gaussianities deserves further 
studies [Tl]. 

Before closing this section we want to comment on an assumption we have implicitly 
made so far. We have imposed on the effective Lagrangian for tt, eq. (jl.3j) . a discrete 
symmetry tt — ► —tt, t — > —t. Without this assumption we should consider an additional 
trilinear operator with the same scaling dimension as tt(Vtt) 2 : 

^(Vtt) 2 V 2 tt, (3.13) 

where 7 is an additional order one coefficient. It is straightforward to repeat the analysis 
in this case; the qualitative behavior of the 3-point function remains the same. 

The tt — * —tt, t —* —t symmetry is equivalent to (ft — > —(ft for the ghost action. In 
absence of this symmetry we can write CPT violating operators in the effective Lagrangian 
[7]. This leads us to a possible interesting effect in the CMBR which is independent of 
how density perturbations are generated: the rotation of the polarization vector of photons 
traveling from the surface of last scattering to us. The origin of this effect is that the ghost 
background {(/)) ^ breaks the CPT symmetry [7j besides the Lorentz one, so that a term 
e^^PF^yAp for the photon field is allowed. This term gives a different dispersion relation 
for the two circular polarizations and it thus rotates the axis of a linearly polarized photon. 
Some limits on the size of this operator come from the observation of distant quasars and 
radio galaxies ^Sj: the rotation must be smaller than roughly one degree in one Hubble 
time. Better limits should come in the future from the CMBR experiments observing 
polarization Temperature fluctuations on the last scattering surface always give a 

certain degree of polarization to the radiation, with a typical "gradient" pattern called E- 
mode. The rotation of the polarization transforms E-modes into "curl" or B-modes giving 
a correlation between temperature fluctuations and B-modes which would otherwise be 
zero. While the present sensitivity is comparable to the limits from distant quasars and 
radio galaxies, planned experiments could reach a sensitivity of 10~ 3 degrees in an Hubble 
time [TU]. We stress that the size of the CPT violating operator cannot be predicted as it 
depends on the coupling of the ghost sector to ordinary matter. 

4 Tilting the potential 

So far we have assumed that the inflationary potential for (ft is exactly flat, i.e. that the 
shift symmetry for (ft is not broken except at (ft = 0. As we have emphasized, unlike slow-roll 
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inflation, no tilt in the potential is needed either to end inflation or to generate acceptable 
density perturbations. However, given that the shift symmetry on cj) must be broken in 
any case, there may well be a tilt in the potential for <j> < 0, although it is technically 
natural for this tilt to be as small as we like. It is thus interesting to see how our previous 
results are modified in the presence of a non-zero slope. We expect, in the limit of very big 
tilt, to recover the predictions of standard slow-roll inflation, where the classical motion of 
4> is dominated by the potential term. As a first approximation, we will assume that the 
potential has a constant slope V = const. The equation of motion for the homogeneous 
perturbations ir then becomes 

7f + 3H-fr + V' = (4.1) 
and we quickly reach the solution for which 

V 

* = -Sff" m 

We see that, because of the potential, we have changed the velocity of the field <f>\ for 
V < 0, we have increased the velocity 1 . Note that, in order to stay within the regime of 
validity of the effective theory, the velocity of tt should still be small compared with M 2 ; 
this motivates us to define the small parameter 

-V 

5 2 = —= . (4.3) 

3HM 2 V ; 

Now, with this non-zero value of tt in the background, the dispersion relation for tt is 

modified by the cubic interaction 

(Vtt^ 2 



so we have (we take a = (3 = 1) 



As long as the k 4 term dominates down to frequencies u ~ H, the 5 term can be justifiably 
neglected; this happens as long as 

5<{H/Mfl 2 . (4.6) 

For larger 5 the new k 2 term dominates before the modes freeze out, at a cross-over 
frequency a> cross determined by 

uw ~ 5 2 M . (4.7) 

We can estimate the size of the perturbations in this case easily, as it is now just like the 
usual case for standard inflation, with the only difference that the spatial momenta are 
rescaled as k — > 5 ■ k. Thus, 

x We will not study the other possibility, V > 0, which leads to a wrong sign spatial kinetic term for -k 

□ 
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which reduces to the usual result of standard inflation 5ith ~ H for S — ► 1, as we expected. 
For the density perturbations we get 

(49) 

p 4> M 2 5 3 / 2 ' K ' ' 

So, ranging over all values of 5, 5p/p goes from (H/M) 5 / 4 (see eq. fOty) to (H/M) 2 . 

Let us now estimate the corrections to the spectral index. It is straightforward to check 
that corrections coming from the variation of H are completely negligible: the variation of 
the potential in one Hubble time is of order AV ~ 5 2 M 4 , which is very small with respect 
to the total vacuum energy Vq. Bigger corrections can come from the variation of <fi in 
eq. (|4.9j) . The effect is of order V" /H 2 and limits on V" come again from the validity of 
the effective field theory description 

5 2 M 2 H ~ \V'\ > \V"\ ■ Acf) = \V"\ ■ (M 2 /H)N e \V"\ < 5 2 H 2 /N e , (4.10) 

where N e is the typical number of e-folds. So that deviations from the flat spectrum can 
be as big as 5 2 /N e : 

K - 1| < S 2 /N e . (4.11) 

The biggest contribution to the tilt comes however from the variation of 5. Assuming a 
variation of 5 during observable inflation of order of 5 itself, we get a contribution to the 
tilt of order 

K - 1| < 1/N e . (4.12) 

Therefore, in the regime 5 > {H/M) 1 / 2 , a rather big deviation from scale invariance is 
possible. 

We can finally also estimate the dimensionless size of the non-Gaussian effects, which 
still come from the same cubic interaction (|3.1[) as before. For the frequencies across < 
uj < M where the k 4 quadratic spatial kinetic term dominates, this operator still has 
scaling dimension 1/4, but for H < lo < u> cmss , the k 2 quadratic term dominates and the 
operator has its familiar scaling dimension of 2. Therefore, the dimensionless size of the 
non-Gaussianity in this case is 

Across \V4 ( H V H 2 6 P 1 



NG ~ CTOSS x — — z£ x — (4 13) 

V M ) UcrosJ MHV 2 p 5 2 1 6) 

so we see that the non-Gaussianities are always parametrically enhanced relative to Sp/p. 
This is to be contrasted with slow-roll inflation, where the leading non-Gaussianity is in fact 
parametrically suppressed relative to 5p/p by slow-roll factors. The same result eq. (|4.13j) 
can be obtained directly evaluating the importance of the non-linear interaction 7r(V7r) 2 
with respect to free field terms at horizon crossing. 

We note that the result for <5 — ^ 1 is what we get in conventional slow-roll inflation, 
with the addition of higher derivative terms suppressed by a cut-off scale M ~ (j) 1 / 2 |12j . 
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6 < {H/M) 1 ' 2 


{H/M) 1 / 2 <5<l 


5~ 1 


Standard slow-roll 


K - !| 





We 


1/iVe 


l/N e 


NG 


ict 2 


io-^- 2 


10" 5 


10" 5 e 



Table 1: Order of magnitude predictions for the tilt in the scalar spectrum and the level 
of non-gaussianity. First column: ghost inflation limit with flat potential. Second: ghost 
inflation with non-flat potential. Third: threshold region, higher derivative terms are as 
important as the tilt in the potential. Fourth: usual slow-roll inflation, higher derivative 
terms are irrelevant. 

This limit can in fact be approached starting from a conventional slow-roll model with 
the addition of higher-dimension operators compatible with the shift symmetry of the 
inflaton. These additional operators will be suppressed by a typical scale M (note that 
now we do not have <fi ~ M 2 , because the motion of the inflaton is not dominated by 
these higher derivative terms). The slow-roll picture makes sense for M 2 S> 4>, otherwise 
all the operators become relevant for the classical motion of the inflaton. In the limit 
M 2 ~ <j) c± V'/H we lose control of the theory, but we qualitatively approach our ghost 
model for <5 c± 1; higher dimension operators become as important as the tilt in the potential 
to describe the classical motion of the inflaton. 

In general if we consider higher derivative interactions suppressed by a typical scale M 
and a potential tilt V', we can look at our model and conventional slow-roll inflation as the 
two extreme limits M 2 ~ S> V'/H (the potential tilt is irrelevant) and M 2 3> 4> — V'/H 
(higher derivative terms are irrelevant). In both the limits we have a trustworthy effective 
field theory description and they merge for V'/H ~ M 2 , i.e. 6 ~ 1, when the potential 
tilt and the higher derivative terms give a comparable contribution to the inflaton motion. 
In this regime both descriptions break down. In the table we summarize the predictions 
for the tilt in the scalar spectrum and the level of non-gaussianity going from the ghost 
inflation limit to the usual slow-roll case. 

5 Conclusions 

We have presented a new way of generating an inflationary de Sitter phase using ghost 
condensation. The spectrum of density perturbations is sharply predicted to have n s = 1 
up to unobservable corrections; gravitational waves are completely unobservable. The non- 
Gaussianities are expected to be large enough to be observed. Unfortunately because there 
are two relevant parameters in the ir effective Lagrangian we can not completely fix the 
magnitude of the 3-point function from the observed density perturbations; however, for all 
parameters of 0(1), the size of the effect is at the level of the current WMAP sensitivity. 
Aside from the overall amplitude, the 3-point correlator ^^(^1)^(^:2)^(^3)^ as a function 
of the spatial momenta k\, fo, &3 is completely determined. 
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We conclude discussing some open issues. As we stressed, all the predictions of our 
scenario are independent of the details of the reheating process, analogously to standard 
slow-roll inflation. But interesting processes might occur at the end of ghost inflation. In 
the region where the shift symmetry is broken we expect that a potential for is generated 
at some level. As discussed in section |1J this potential will kick <p out of its stable value 
<j> = M 2 and (j> will be driven back to its original value by the expansion of the universe once 
the shift symmetry is restored. The stress energy tensor of this homogeneous perturbation 
7r redshifts like matter [7j (p oc a~ 3 ) so that it would be interesting to study the evolution 
of this fluid to see if it could be a viable candidate for dark matter. 

Another interesting possibility is that inflation ends with a phase transition, after which 
the ghost is not condensed anymore: Lorentz symmetry is restored. This scenario allows 
to raise the scale M, because we are not limited anymore by the present modification 
of gravity. Obviously the phase transition leads us outside the regime of validity of the 
effective field theory around the ghost condensate, but it is likely that all our results about 
density perturbations still apply. 
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Appendix: conservation of £ 

The use of the formula (|2.1|) for the calculation of the gauge- invariant quantity £ is justified 
in general by the conservation of £ during the reheating process. In this appendix we shall 
explicitly check the conservation. For simplicity we take the limit in which the time scale 
associated with reheating is much shorter than the Hubble time scale and assume that 
spatial components of the stress-energy tensor remain non-singular. In this case we can use 
the Israel junction condition without a singular matter source to obtain the matching 
condition for perturbations even without any knowledge about the details of microscopic 
processes. 

Starting with general scalar-type perturbations in the longitudinal gauge, 
ds 2 = - [1 + 2$(i, a?)] dt 2 + [1 - 29(f, x)] a 2 (t)S ij dx i dx j , <j> = c/> (0) (t) + vr(t, x) , (A.l) 
we can perform an infinitesimal gauge transformation 

Sg^ - Sg^ - e w „ - , vr - 7T - pd^W , (A.2) 
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where 



>°)(*o) 

-a\t) 



so that the metric and the scalar field become 
ds 2 = —dt 2 + qijdx l dxi , 



^{t > ,x)dt > 



t 



4 



to 



a?{t>) 



(A.3) 



(1 - 2* + 2#£ )a 2 <% - diij - djti , 



(j, = ^(t)+ir(t,x)-Tr(t ,x)-^°\t) <S>(t',x)dt' 



(A.4) 



to 



Here, H = a/a and to is the value of t when the homogeneous background 4>(°\t) reaches 
the critical value 4>q at which the phase transition is triggered. The coordinate system after 
the gauge transformation is nothing but the Gaussian normal coordinate system based on 
the reheating hypersurface 4> = 4>q. In the Gaussian normal coordinate system, the Israel 
junction condition ^Tf\ at the reheating hypersurface is 



\X] = lim X(t) - lim X(t) . 
t-vto+O f->-to-0 



Thus, we obtain 







for the unperturbed geometry, and 

; [*] = ; 



7T 



[H}=0 



7T 







(A.5) 
(A.6) 

(A.7) 



(f)(0) 

for the linear perturbation. 

Since in a pure de Sitter background the gravitational potential at large scales decays [7j, 
it becomes vanishingly small after a long period of de Sitter expansion. Hence, we set 
$ = \£ = for t < tQ. On the other hand, the universe after reheating is dominated by 
radiation. In the radiation dominated epoch, it is known that the Fourier components 
and *$>k of $ and ^, respectively, at superhorizon scales behave as 



^ k ~ A k + B k 



a(t ) 



(t > t ) 



(Ai 



where A k and B k are constants. The coefficients A k and B k are determined by the Israel 
junction condition as 

A k = -B k = , (A.9) 

6 r ' t=to 
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where iTk(t) is the Fourier component of ir. Therefore, we obtain the following behavior of 
and ^ after reheating. 



1 



a(t ) 
a(t) 



2 7Th 

-H- ' 



3 0(0) 



(A.10) 



t=t 



Finally, with the help of the well-known relation £ = — 1$ in radiation-dominated epoch, 
this shows that 



7T 



0(0) 



(A.ll) 



and that we can safely use the formula (|2.1f) . In the main part of this paper we shall omit 
the superscript (0) for rf>(°> . 
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